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LONG-TERM  GOALS 

To  develop  a  robust  multi -phase,  multi-class  numerical  modeling  framework  for  both  cohesive  and 
non-cohesive  sediment  transport  in  the  fluvial,  estuarine  and  coastal  environments. 

OBJECTIVES 

This  study  specifically  focuses  on  numerical  modeling  of  critical  processes  at  small-scale  (0(cm- 
100m)).  Specific  objectives  are: 

1 .  To  develop  turbulence-resolving  numerical  model  of  fine  sediment  transport  in  the  oscillatory 
boundary  layer  in  order  to  understand  how  turbulence-sediment  interactions  can  determine 
fluid  mud  transport  and  the  state  of  muddy  seabed. 

2.  To  develop  a  two-dimensional-vertical  (2DV)  numerical  model  based  on  Reynolds-averaged 
Navier-Stokes  (RANS)  equations  and  volume  of  fluid  method  for  free-surface  tracking  to  study 
mechanisms  causing  landward  and  seaward  fine  sediment  transport  in  inter-tidal  zones. 

3.  To  develop  2DV-RANS  modeling  for  wave-mud  interactions  in  order  to  understand  the 
competing  effects  between  mud  dissipation  and  shoaling  in  detennining  the  resulting  nonlinear 
wave  propagation. 

APPROACH 

Cohesive  sediment  transport  involves  a  variety  of  physical  mechanisms  including  boundary  layer 
processes  (tidal  and  wave),  gravity-driven  flow,  turbulence  modulation,  flocculation,  non-Newtonian 
rheological  behavior  and  consolidation  (e.g.,  Mehta  1989;  Winterwerp  and  van  Kerstem  2004).  A 
general  modeling  framework  appropriate  for  a  wide  range  of  concentration  needs  to  be  based  on 
multiphase  flow  theory.  In  this  study,  a  fine  sediment  transprt  modeling  framework  based  on 
Equilibrium  Eulerian  Approximation  (Ferry  &  Balachandar  2001)  to  the  multiphase  equations  has  been 
developed  and  extended  to  model  various  cohesive  sediment  transport  processes.  This  fine  sediment 
modeling  framework  is  the  basis  of  three  numerical  models  developed  for  1DV  (one-dimensional- 
vertical),  2DV  and  3D  for  different  applications  investigated  in  this  study.  In  the  Reynolds-averaged 
1DV  modeling,  the  dynamics  of  wave-supported  gravity-driven  mudflows  has  been  studied  (Hsu  et  al. 
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2007;  2008).  Efforts  are  also  devoted  to  improve  existing  flocculation  formulation  with  more  robust 
parameterizations  on  fractal  dimension  and  floe  yield  strength  (Son  &  Hsu  et  al.  2008,  2009).  New 
flocculation  formulation  is  then  adopted  in  the  1DV  model  to  study  the  effect  of  flocculation  and  bed 
erodibility  on  cohesive  sediment  transport  in  a  meso-tidal  environment  (Son  &  Hsu  20 10a, b).  In  the 
Reynolds-averaged  2DV  modeling,  we  further  utilize  a  volume  of  fluid  (VOF)  code  for  free-surface 
tracking  (Lin  &  Liu  1998)  and  investigate  wave  propagation  over  muddy  seabed  (Torres-Freyermuth 
&  Hsu  2010).  Recently,  this  2D V- VOF  code  is  also  extended  to  simulate  cross-shore  fine  sediment 
transport  in  intertidal  flats.  To  greatly  reduce  the  uncertainties  in  the  turbulence  closure  for  sediment¬ 
laden  flow,  a  3D  highly  accurate  turbulence-resolving  Navier-Stokes  solver  is  extended  to  study  the 
role  of  sediment-induced  density  stratification  in  determing  the  resulting  fine  sediment  transport  and 
bed  state  in  the  wave  boundary  layer  (Ozdemir  et  al.  2010a,  b). 

WORK  COMPLETED 

1DV  modeling  of  wave  supported  gravity-driven  mudflow  and  the  effects  of  flocculation  on  fine 
sediment  transport  in  tide-dominant  environment. 

These  parts  of  the  investigation  have  been  completed.  Main  findings  are  summarized  in  the  report  of 
FY09.  Five  journal  papers  are  published,  in  press  or  currently  under  review. 

2DV  RANS-VOF  modeling  of  wave-mud  interaction 

Ongoing  studies  supported  by  ONR  to  study  wave-mud  interaction  at  Louisiana  shelf  have  revealed  the 
role  of  direct  dissipation  and  nonlinear  wave-wave  interactions  in  determining  observed  wave 
attenuation  (Sheremet  &  Stone  2003;  Elgar  &  Raubenheimer  2008).  Our  ongoing  modeling  work  using 
2DV-RANS-VOF  wave  model  suggests,  on  horizontal  bottom,  low  frequency  wave  energy  similar  to 
that  of  high  frequency  waves  (Sheremet  &  Stone  2003),  is  transferred  to  sea-swell  band  and  get 
dissipated  (Torres-Freyermuth  &  Hsu  2010).  On  the  other  hand,  for  the  case  of  sloping  seabed, 
shoaling  effect  dominates  and  drains  energy  from  sea-swell  to  low  frequency  band,  consistent  with 
field  observation  of  Elgar  &  Raubenheimer  (2008).  Our  model  results  suggest,  however,  comparing 
with  results  without  mud,  the  presence  of  mud  can  still  reduces  the  amount  of  energy  transferring  to 
low  frequency  band  due  to  shoaling.  Hence,  there  exists  a  competition  between  direct  mud  dissipation 
and  shoaling  process  at  muddy  shelf.  This  problem  is  not  only  dependent  on  mud  rheology  but  also 
flow  depth  and  beach  slope  that  deserves  further  numerical  investigations.  Essentially,  without  a 
complete  understanding  on  the  competing  effects  between  bottom  mud  dissipation  and  shoaling  in 
determining  the  wave  energy  evolution  and  dissipation,  it  is  difficult  to  back  calculate  the  bottom 
seabed  properties  from  the  observed  wave  field.  A  journal  publication  summarizing  our  numerical 
model  development  effort  is  published  (Torres-Freyermuth  and  Hsu  2010).  More  recent  investigation 
on  the  competing  effects  between  bottom  mud  dissipation  and  shoaling  over  sloping  seabed  is 
underway.  Numerical  investigations  are  carried  out  for  different  beach  slopes  and  mud  availability. 

3D  two-way  coupled  turbulence-resolving  simulation 

A  3D  numerical  simulation  tool  for  fine  sediment  transport  in  oscillatory  boundary  layer  has  recently 
been  developed  where  all  the  scales  of  carrier  flow  turbulence  are  resolved  without  the  uncertainties  in 
turbulence  closure  (Ozdemir,  Hsu  and  Balachandar,  2010,  J.  Fluid  Mechanics,  in  press).  For  a  given 
Stokes  Reynolds  number  and  settling  velocity,  we  vary  the  available  amount  of  sediment  in  the 
oscillatory  boundary  layer  and  observe  four  different  flow  regimes  ranging  from  well-mixed  sediment 
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in  the  wave  boundary  layer,  formulation  of  lutocline  and  laminarization.  Our  ongoing  work  focuses  on 
investigating  the  existence  of  these  flow  regimes  and  the  lutocline  dynamics  for  different  settling 
velocities  and  Reynolds  number. 

2DV  RANS-VOF  modeling  of  cross-shore  sediment  transport  at  tidal  flats 

Conventional  coastal  modeling  systems  have  difficulties  in  resolving  the  tidal  front  at  the  wetting  and 
drying  seabed  and  some  numerical  approximations,  such  as  specifying  a  minimum  artificial  flow 
depth,  are  often  adopted.  However,  it  is  not  clear  how  such  artificial  numerical  treatment  can  affect  the 
accuracy  of  the  results  when  tidal  water’s  edge  is  of  interest.  Hence,  a  new  numerical  modeling 
approach  for  this  problem  is  necessary.  The  2DV  RANS-VOF  numerical  model  previously  used  for 
wave-mud  interaction  is  revised  for  tidal  flow.  Our  preliminary  model  results  suggest  the  VOF 
approach  is  able  to  better  predict  the  wetting  and  drying  processes  and  resolve  the  turbid  tidal  water’s 
edge.  The  only  empirical  parameter  in  the  VOF  approach  that  is  sensitive  to  the  water’s  edge  process  is 
the  bed  roughness  height.  However,  bed  roughness  height  is  a  much  more  physically-based  parameter 
that  can  be  measured  with  reasonable  confidence. 

RESULTS 

The  highlights  of  our  study  on  cohesive  sediment  transport  for  FY10  is  summarized  here: 

The  state  of  muddy  seabed  in  wave-dominant  environment 

One  of  the  most  important  findings  resulted  from  AMASSEDS  indicates  a  significant  reduction  of 
bottom  drag  coefficient  as  tidal  currents  propagate  over  muddy  seabed  (Beardsley  et  al.  1995)  due  to 
attenuation  of  turbulence  via  fluid-mud-induced  density  stratification  (Trowbridge  &  Kineke  1994). 

On  the  other  hand,  field  and  analytical  studies  on  wave  propagation  over  muddy  seabed  (e.g.,  Sheremet 
&  Stone  2003;  Dalrymple  &  Liu  1979)  report  increasing  energy  dissipation  due  to  the  presence  of 
bottom  fluid  mud.  The  completely  contradictory  outcome  of  mud  on  the  hydrodynamic  dissipation  for 
(tidal)  currents  and  waves  deserves  more  detailed  investigation.  Through  ongoing  ONR  supports,  we 
have  developed  a  3D  turbulence-resolving  simulation  tool  based  on  highly  accurate  pseudo-spectral 
scheme  to  study  fine  sediment  transport  in  oscillatory  boundary  layer  (Ozdemir  et  al.  2010a,  b).  Our 
simulation  results  reveal  four  distinct  regimes  of  wave-induced  fine  sediment  transport  as  sediment 
concentration  increases  (see  Figure  1;  Simulations  are  of  Stokes  Reynolds  number  Re 5=  1 000  and 
settling  velocity  0.5  mm/s).  These  different  flow  regimes  are  caused  by  sediment-induced  density 
stratification  on  turbulence  and  mixing  in  the  wave  boundary  layer.  Specifically,  we  find  that  in  the 
regime  II  of  moderate  concentration  (10~50g/l),  pronounced  lutocline  features  appear  at  the  top  of  the 
wave  boundary  layer.  However,  no  significant  damping  of  carrier  flow  turbulence  (hence  drag 
reduction)  within  the  wave  boundary  layer  is  observed  because  damping  mostly  occurs  around  the  top 
of  wave  boundary  layer  where  turbulence  production  is  already  small  (i.e.,  density  stratification  merely 
damps  the  nonlocal  transport  terms).  Significant  damping  of  carrier  flow  turbulence  and  laminarization 
are  observed  in  regime  III  of  higher  sediment  concentration  (»50  g/1).  However,  under  such  high 
concentration,  rheological  stress  may  become  important  to  attenuate  wave  energy.  Better 
understanding  on  the  existence  and  main  features  of  these  flow  regimes  are  critical  to  our  predictive 
skill  on  the  bed  state  and  seabed  properties. 
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Mechanisms  causing  landward  and  seaward  sediment  transport  in  a  tidal  flat 

To  understand  the  morphodynamics  of  intertidal  zones,  it  is  critical  to  study  mechanisms  causing 
landward  and  seaward  sediment  fluxes  subjected  to  tidal  flows  and  waves.  The  1DV  model  for 
cohesive  sediment  transport  is  developed  to  quantify  how  sediment  flux  is  associated  with  tidal 
skewness  (asymmetry)  (Son  and  Hsu  2010).  In  the  literature,  the  settling-lag  effect  is  considered  to  be 
the  main  mechanism  causing  landward  transport  (e.g.,  de  Swart  &  Zimmerman  2009).  On  the  other 
hand,  the  export  of  fine  sediment  through  channel-runnel  systems  as  well  as  the  wave  effects  has  been 
identified  as  one  of  the  main  offshore  transport  mechanisms  in  tidal  flats.  Field  observations  in  the 
ongoing  Tidal  Flat  DRI  measure  sediment  pulse  during  flood  and  ebb  tidal  water’s  edge  passage  (A. 
Ogston,  personal  communication),  consistent  with  several  earlier  studies  in  the  European  tidal  flat 
(e.g.,  Christie  et  al.  1999).  This  provides  important  evidence  to  better  understand  landward  and 
seaward  transports  in  details,  especially  in  terms  of  the  contrast  between  the  channel  and  the  flat. 
Estuarine  modelers  in  the  ongoing  Tidal  Flat  DRI  (e.g.,  Chen  et  al.  2010)  also  identify  export  through 
channel  as  one  of  the  missing  mechanism  in  the  existing  coastal  models,  which  prevent  coastal  models 
to  capture  the  short  residence  time  of  river-delivered  fine  sediment.  Because  major  transport  occurs  at 
tidal  water’s  edge  (e.g.,  Christie  and  Dyer  1998),  there  is  a  need  to  further  investigate  this  issue  with 
detailed  numerical  model  that  can  resolve  processes  in  the  very  shallow  depth. 

In  the  present  two-dimensional-vertical  (2DV)  modeling,  we  primarily  focus  on  idealized  study  in  the 
cross-shore  direction  (the  direction  perpendicular  to  the  tidal  water’s  edge).  Flow  and  sediment 
conditions  investigated  here  are  of  a  flat  slope  0.0013,  tidal  range  4  m  (similar  to  Willapa  site), 
sediment  settling  velocity  0.5  mm/s,  critical  bed  stress  0.15  Pa.  Idealized  standing- wave  type  tidal 
forcing  of  sine  shape  is  specified  (tidal  asymmetry  is  zero).  Two  snapshots  of  the  model  results  during 
flood  and  ebb  are  presented  in  Figure  2.  Toward  the  middle  of  the  flood  (left  panel  in  Figure  2),  tidal 
velocity  is  strong  and  significant  amount  of  sediment  is  suspended.  The  tidal  water’s  edge  induces 
strong  velocity  and  sediment  resuspension.  This  turbid  tidal  edge  is  a  commonly  observed  feature  in 
the  intertidal  zones  (Christie  and  Dyer  1998).  Right  after  the  peak  ebb  velocity  (right  panel  in  Figure 
2),  ebb  turbid  tidal  edge  can  be  clearly  seen.  However,  numerical  model  predicts  different  ebb  tidal 
water’s  edge  from  that  during  flood.  The  ebbing  flow  leaves  a  thin  water  layer  tail,  which  makes  a 
more  gradual  transition  of  velocity  at  ebb  water’s  edge.  The  asymmetry  during  the  flood  and  ebb  tidal 
edge  can  be  more  clearly  identified  based  on  the  time  series  of  velocity  profiles  and  sediment 
concentration  profiles  at  x=3900  m  (Figure  3).  The  flood  velocity  induces  a  sediment  concentration 
spike  that  is  of  high  magnitude  but  short  duration.  One  the  other  hand,  the  ebb  velocity  induces  smaller 
concentration  but  with  longer  duration.  Such  asymmetry  varies  in  the  cross-shore  direction  and  gives 
seaward  transport  in  the  lower  flat  and  landward  transport  in  the  upper  flat.  More  detailed  investigation 
is  underway  to  compare  the  present  model  results  with  the  classic  equilibrium  and  settling-lag 
paradigms  (e.g.,  Friedrichs  and  Aubruy  1996). 

IMPACT/APPLICATIONS 

The  present  research  efforts  produce  a  numerical  modeling  framework  for  cohesive  sediment  transport 
for  various  applications  in  tide-  or  wave-dominated  environment.  Our  model  development  efforts  shall 
contribute  the  modeling  capability  in  the  ongoing  ONR  related  research  effort  on  Tidal  Flat  DRI, 
Wave-mud  Interaction,  Community  Sediment  Transport  Modeling  System  (NOPP-CSTMS),  Inlet  and 
River  Mouth  DRI. 
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Figure  1:  Simulation  results  for  dilute  (regime  I),  medium  (regime  II)  to  high  concentration  (regime 
III)  cases  (near  bed  concentration  C~O(0.1),  10,  and  50  g/l,  respectively).  1st  row:  time-series  of 
free-stream  velocity  and  associated  phases;  The  2nd  and  3rd  rows:  iso-surface  of  sediment 
concentration  for  regime  II  (C~10  g/l),  and  regime  III  (C~50  g/l)  at  three  different  phases.  When 
concentration  is  large,  flow  laminarizes  under  wave  crest  but  shear  instability  occurs  during  flow 
reversal  (see  3rd  row).  These  features  are  very  different  from  that  of  lower  concentration  where  flow 
remains  turbulent  throughout  the  wave  cycle  (see  2nd  row).  The  4th,  5th,  and  6th  rows  are  averaged 
concentration,  streamwise  velocity  and  RMS  turbulent  velocity  fluctuations  profiles.  Significant 
damping  of  turbulence  due  to  sediment-induced  density  stratification  causes  laminarization  of 
velocity  and  reduction  of  wave  boundary  layer  thickness  (compare  blue  with  other  curves).  Magenta 
symbols  in  the  5th  row  are  theoretical  laminar  solution.  In  regime  III,  flow  is  almost  laminarized  by 
mud  except  during  flow  reversal.  When  concentration  is  0(100)  g/l  or  more  (Regime  IV),  flow  is 
completely  laminarized  (not  shown  here).  Re  §=1000;  Ws/U=0.0009. 
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x  (m)  x  (m) 

Figure  2:  Snapshots  of  numerical  model  results  during  (a)  flood  tide  and  (b)  ebb  tide  for  flow  field 
and  sediment  concentration  (upper  panel)  and  cross-shore  distribution  of  flow  velocity  (dashed-blue 
curves)  and  sediment  concentration  (solid-red  curves)  at  2.5  cm  above  the 

bed  (lower  panel). 
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Figure  3:  Time  series  of  vertical  profiles  of  cross-shore  velocity  magnitude  (upper  panel),  bottom 
stress  (middle panel)  and  sediment  concentration  (lower  panel)  at  x=3900  m.  The  red-dashed  line  in 
the  middle  panel  represents  the  critical  bed  stress  (0.15  Pa)  specified 

in  this  case. 
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